NACA TM 1420 


■ - 


-4 

CD 

->4 

O 


r^? o 
* * / 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL MEMORANDUM 1420 


ON THE CONTRIBUTION OF TURBULENT BOUNDARY 
LAYERS TO THE NOISE INSIDE A FUSELAGE 

by 

G. M. Corcos 
and 

H. W. Liepmann 
Douglas Aircraft Company, Inc. 



Washington 
December 1956 


A rsf rvo 

t: ; “ '■ ... u 


*■ t - « 
i . . . 


-1. S " 1 / t 1 


□ 

fcH 

-C 

-C 

03 

UJ 


- -'V 

!. I 


I 


TECH LIBRARY KAFB, NM 


TECHNICAL MEMORANDUM lk20 


TECH LIBRARY KAFB, NM 



i mil i 

01444B3. 


ON THE CONTRIBUTION OP TURBULENT BOUNDARY 
LAYERS TO THE NOISE INSIDE A FUSELAGE* 


by 

G. M. Coreas'** 
and 

H* W* Liepmann*** 


September, 1956 


"Unedited by the NACA (the Committee 
takes no responsibility for the 
correctness of the author's 
statements . ) M 

x % 

Aerodynamics Research Group 
Douglas Aircraft Company, Inc, 

Santa Monica, California 

* X California Institute of Technology 
Pasadena, California 
(Consultant, Aerodynamics Research 
Group, Douglas Aircraft Company, Inc. 
Santa Monica, California) 




NA.CA TM 1420 

TABLE OF CONTENTS 

m 

. ABSTRACT I 

INTRODUCTION 1 

I. THE ACOUSTIC COUPLING OF A RANDOMLY 

VIBRATING PLATE WITH AIR AT REST k 

II. THE DYNAMIC BEHAVIOR OF THE SKIN 8 

a) The Mean Acceleration 9 

h) The Length Scale A 12 

III. THE FORCING FUNCTION ik 

IV. SPECIAL CASES 20 

1. Convected Turbulence 20 

a) The coupling of the plate -with air at 

rest in the case of convected turbulence • • . . 20 

b) The response of the plate 21 

2. The case of zero scale 23 

V. SUMMARY OF RESULTS AND DISCUSSION 2k 

A Note on Testing 25 

VI. REFERENCES 2TJ 

APPENDIX I: THE RANDOM RADIATION OF A PLANE SURFACE: 

A. Four Limiting Cases 28 

B. The Noise Generated by Skin Ripples 

of Fixed Velocity 35 

C. The Generation of Noise by Plate 

Deflections of Zero Scale 38 

APPENDIX II: THE SIMPLIFICATION OF THE PLATE RESPONSE INTEGRAL . . kO 

■t 

APPENDIX III: THE EVALUATION OF THE INTEGRAL SCALE 42 




NA.CA TM 1420 


1 


ABSTRACT 


The following report deals in preliminary fashion with the transmission 
through a fuselage of random noise generated on the fuselage skin "by a tur- 
bulent boundary layer. The concept of attenuation is abandoned and instead 
the problem is formulated as a sequence of two linear couplings: the turbu- 

lent boundary layer fluctuations excite the fuselage skin in lateral vibra- 
tions and the skin vibrations induce sound inside the fuselage. The techniques 
used are those required to determine the response of linear systems to random 
forcing functions of several variables. A certain degree of idealization has 
been resorted to. Thus the boundary layer is assumed locally homogeneous, 
the fuselage skin is assumed flat, unlined and free from axial loads and the 
"cabin" air is bounded only by the vibrating plate so that only outgoing 
waves are considered. Some of the details of the statistical description 
have been simplified in order to reveal the basic features of the problem. 

The results, strictly applicable only to the limiting case of thin 
boundary layers, show that the sound pressure intensity is proportional to 
the square of the free stream density, the square of cabin air density 
and inversely proportional to the first power of the damping constant and 
to the second power of the plate density. The dependence on free stream 
velocity and boundary layer thickness cannot be given in general without a 
detailed knowledge of the characteristics of the pressure fluctuations in 
the boundary layer (in particular the frequency spectrum). For a flat 
spectrum the noise intensity depends on the fifth power of the velocity 
and the first power of the boundary layer thickness. This suggests that 
boundary layer removal is probably not an economical means of decreasing 
cabin noise. 

In general, the analysis presented here only reduces the determination 
of cabin noise intensity to the measurement of the effect of any one of 
four variables (free stream velocity, boundary layer thickness, plate 
thickness or the characteristic velocity of propagation in the plate). 

The plate generates noise by vibrating In resonance over a wide 
range of frequencies and increasing the damping constant is consequently 
an effective method of decreasing noise generation. 

One of the main features of the results is that the relevent quantities 
upon which noise intensity depends are non-dimensional numbers In which boundary 
layer and plate properties enter as ratios. This is taken as an Indication 
that in testing models of structures for boundary layer noise It is not 
sufficient to duplicate In the model the structural characteristics of the 
fuselage. One must match properly the characteristics of the exciting 
pressure fluctuations to that of the structure. 
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INTRODUCTION 


In his efforts to minimize the noise levels for which he is responsible , 
the airplane designer has had to pay increasing attention to a source of 
noise which until recently had been ignored. This is the boundary layer. 

The boundary layer will generate noise whenever it is the seat of any 
fluctuating phenomenon. In particular it will nurture random pressure 
fluctuations whenever it is turbulent. 

The designer* s interest will naturally center on the characteristics 
of that part of the boundary layer noise which has been transmitted through 
the fuselage skin and into the cabin rather than on that part which is 
radiated into the free -stream. This is so because the radiation intensity 
is, as we shall see, a rapidly increasing function of the velocity of the 
boundary with respect to the air and to a likely observer outside a fuselage 
either the relative velocity of the plane is low (as near a take-off or 
landing) or the plane is considerably distant. 

As a consequence the practical question which has to be raised concerns 
the effects on a fuselage skin, and on the air which it encloses, of the 
boundary layer pressure fluctuations acting directly on the skin.* 

Two features of this problem are worth noting: To begin with, the 

fuselage will transmit noise only by deflecting laterally. The thickness 
of the skin is very small when compared to the wave length in metal of 
audible sound waves so that, effectively there are no fluctuating pressure 
gradients within the skin and hence the latter will not oscillate in lateral 
compression. 

In the second place, the turbulent pressure fluctuations in the boundary 
layer are random both in space and in time. The fluctuations are generated 
locally. If they are measured simultaneously at two different points of 
the boundary layer, say on the skin itself, they are found to have no 
relationship with each other unless the two points are separated by a very 
short distance. Alike the value of the pressure fluctuation at a given 
point soon loses correlation with itself.** 


*Noise intensity is defined in this report as PI 2 /p 0 a i* It is assumed that 
this is the physical quantity of interest. It has the dimensions of an 
energy flux; but it is not necessarily equal to the energy flux at some 
point in the field, nor is it necessarily equal to the density of 
energy radiated away (lost) by the fuselage. 

^Recently two authors (Refs. 5 & 6) have suggested that the randomness in 
time is not independent of the randomness in space; i.e., that the pressure 
fluctuations at the wall are created by the convection at a single speed 
of a "frozen" pattern of pressure disturbances. Some attention is paid 
later to this eventuality which is treated as a special case. 
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We are thus led to visualize the process of transmission of boundary- 
layer noise through the metal skin of the fuselage as follows: A multitude 

of external pressure pulses push the elastic skin in and out and the skin, 
in turn, not unlike a set of distributed pistons, creates inside the fuse- 
lage pressure waves which propagate and superimpose. This constitutes cabin 
noise . 

It is of course desirable to determine the characteristics of this 
noise. One should point out that the data for the problem are not complete 
and are not likely to be so in the near future. Specifically the structure 
of the turbulent mechanism within the boundary layer and, in particular 
of the coupling between pressure fluctuations, velocity fluctuations and 
temperature fluctuations is not enough explained or measured to define 
wholly our forcing function. As a consequence it is not now possible to 
define say average cabin noise intensity p-^ 2 as a function of say, free 
stream Mach number, Reynolds number and plate characteristics. It Is 
however possible and it is the purpo se of this report to indicate the 
approximate functional dependence of p^ on these quantities and thus 
to give similarity rules which will reduce to a minimum the amount of 
testing required. 

We assume at the outset that the boundary layer unsteady pressure 
field is known and that it induces small deflections In the skin. As a 
consequence 

(a) The skin dynamics are described by a linear equation 

(b) The generation of a random pressure field inside the cabin 
is a linear radiation problem. 

Thus the mathematical techniques used are those required to obtain 
the response of linear systems to stochastic forcing functions of several 
variables . 

We also assume the fuselage to be a large flat plate. This assumption 
is not necessary but it simplifies the discussion and allows us to present 
more clearly the new features of the problem. 

The material in this report is presented as follows: First we study 

the radiation of sound from a randomly vibrating plate. It Is found that 
the sound levels in the cabin are defined by the intensity and the scales 
of the plate normal accelerations. 

Second, generalized Fourier analysis is put to use in order to relate 
the normal acceleration of the plate to the forces exerted "on it by the 
boundary layer flow. 

Third, the boundary layer forces are defined in terms of flow character- 
istics, and dimensional similarity is used to determine the significant 
parameters . 

Finally, the functional form of the noise intensity in the cabin is 
given save for an unknown function of one non-dimensional parameter. This 
function depends on the frequency spectrum of the pressure fluctuations 
in the boundary layer. Wo measurements yielding this spectrum have been 
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reported to date and speculations concerning it would introduce in the 
analysis both complication and uncertainty. 

A summary of results is given at the end of the report. Some 
derivations and some of the longer arguments have been presented as 
appendices to the text. 
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I. THE ACOUSTIC COUPLING OF A RANDOMLY 
VIBRATING PLATE WITH AIR AT REST 


We start with Rayleigh's well known solution of the acoustic equations 
when the sound is generated in an otherwise “unbounded stationary gas by a 
large flat plate or disc oscillating normally to its plane (Ref. 1 page 107). 



where : 



(i) 


the static pressure p has been broken into a steady part p and 
a fluctuating part , p^ : 

P = P* (x,y,z) + 2 ± (x,y,z,t) 

V n = the normal velocity of the plate 

do* = an element of plate surface area 

a^ = the speed of sound in the air 

3 = the vector position over which the integration is carried out 
= (x',o,z') 

1*1 = the distance between source and field points 
= \f~ (x-x') 2 + y 2 + (z-z*) 2 = r 

A = the total area of the plate 


The subscript i refers to properties inside the cabin or on the side of the 
plate on which air is not flowing. 

The normal acceleration ^Vn/Jt is a random function of time and space 
and so is p^. We wish to evaluate the mean square of the pressure (a quantity 
to which sound intensity is directly proportional) . 


We notice that for a sum 

( 'Z <2i) y 

* 2 at f aj 

alike, for the integral in (l)* 




*This discussion follows closely the arguments set forth in ref . 3 and the 
background of information on statistical methods can be found in ref. 2. 
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^ ' ( I J„ W(^ t- i) i A I f ( ^' c - ’S) ) 


and since ax = ax if a does not depend on time , we can write 


>L 


r 

4tc 2 J 

Aft 

Nov the expression. 




dt 


( 2 ) 


a: ' a, 1, 



is the correlation (i.e. the time average of the_productJ> of the normal 
acceleration of the plate surface at two points and So and at two 
different times (t- ) and (t - If the two points vibrate 

completely independently from each other, this expression is zero and 
if the two points are brought together ( Sj -»■ So * therefore the 

correlation function is simply (dVn/ a f) z - Now we 'assume that the average 
properties of the plate motion are the same anywhere on its surface and 
at any time (we assume statistical homogeneity and stationarity) - Then the 
expression above for the correlation becomes merely a function of the 
distance ( \ S^- Sgl) between the two points and of We call this 

function m : 


then 






^‘ t ~. i d. vi~ j. 


( 3 ) 


Now we can evaluate pj_2 (Y) under a variety of assumptions for7/j , for the 
plate area A end for the distance Y between the plate and the observer. We 
will consistently hold the view that the normal acceleration at most points 
of the plate surface are not correlated O) and that two points of 
the plate, in order to show appreciable correlation, must be a small 
fraction of the total plate size away from each other. 


We define as A- the mean distance over which i^yk) is strongly 
correlated (i.e. 3M[jggM) and call it the integral (length) scale* , 

* The integral scale is given, say in the x-direction by 

( *«•/*)* » r (v^o w dxz 
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Our hypothesis can then he expressed as 

K « H 

where R is the average linear dime nsion of the plate. 


A representative case 

Suppose that the observer distance Y to the plate is such that 

\ « y « 

and that no appreciable phase difference can arise at Y between two strongly 
correlated signals (i.e., between sound pulses originating within \ of 
each other). Then, if we examined equation ( 3 ) 


$( 4£j oL<J - v '! d.«-> 

- 1 

we see tjoat the inner integral contributes very littl| except when the 
point b t is appr oximate ly within a distance X. of S, • Then, approximately 
r\^m > 1 | and Y s iSYnTTR Thus the inner integral is approximately 

^ t dt ) 



and equation ( 3 ) can be rewritten 



We should note that 

The integral in (4) is not finite if the plate area A is infinite 

The length scale \ plays no role in the geometry of the problem. 

The integral is a function of Y, the distance to the plate, and the plate 
dimensions only. For instance if the plate is circular and of radius R 


* This case is treated more rigorously in Appendix IA 
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U -i 
r 


It is apparent from this result that the distance Y from the plate to 
the observer is measured in terms of plate diameters, and not in terms 
of average correlation lengths or wave lengths \ . This result holds for 
all the cases considered (see Appendix I) and depends only on the 
assumption that at a given time the plate vibrations are largely un- 
correlated or incoherent.* 


Here, in effect, 
to define a strength. 


loses Its Identity and combines with 




Equation (^) indicates that in order to evaluate the pressure ■ 
intensity for the case considered above, one must first determine (oVfy0£) 
the mean square normal acceleration of the plate and ?\ , the length 
scale for the plate deflections. Other cases (i.e. cases for which 
either the plate deflects differently or the observer moves closer to 
it) are treated in Appendix IA. For seme of them the time scale or 
mean period 10* is required as well as \ . 


We rewrite Equation (4) 






( 6 ) 


where C 9 ' V ) is a function of the plate geometry and of the distance V 
between the observer and the plate. It Is defined from equation (4) as 

- Ac* 

Notice that if the integral scale is not the same in the x* and in 
the z* direction we may simply substitute in equation (6) and AV 

for X 


* This result holds, as Appendix IB shows, even when the pressure is 

generated by the travel through air at rest of a "randomly bumpy plate" 
i.e., when the time-wise and space-wise variations of upwash are not 
Independent. The latter example is therefore quite distinct from the 
flow of an Infinite (periodically) wavy wall for which the only 
characteristic length is the wave length. 
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II. THE DYNAMIC BEHAVIOR OF THE SKIN 


The response of the bare fuselage skin to the random pressure field 
of the boundary layer is given, in the absence of axial loads, by a plate 
equation. For a flat plate this equation is 


Eh* 

ia-O-/-') 



Itt * fit ‘i 


(7) 


where x and z are the coordinates along the plate surface (x being the 
free stream direction), y the deflection of the plate at a point, E the 
modulus of elasticity of the plate material, o- its density, Poisson's 
ratio, 2h the thickness of the plate, /3 a damping constant which has 
dimensions (l/T). Damping may be present because energy is absorbed either 
within the skin or by the air. For air damping 

Notice that to a, How for air damping is to provide for a feedback in 
the coupling between the plate and the air at rest. On the other hand 
we exclude feedback between the plate and the boundary layer. In other 
terms we are not considering the possibility that the plate vibrations 
are large enough to induce time -dependent pressure gradients of the same 
order of magnitude as our forcing function. Such a feedback would amount 
to panel flutter. It cannot be handled by the present method. 


f(x,z,t) is the random force / unit mass exerted by the pressure 
fluctuations on the plate surface. It is characterized by a power spectral 
density F(k 1 ,k2jt J ) which is a continuous function of the wave numbers kj_ 

(in the x direction), kg (in the z direction) and of the. frequency CD. 

The coefficient E/3<r(l-/^) has the dimensions of a velocity squared 
and it is defined as c^. 


*The skin construction may be such that the damping it causes is primarily 
viscous or primarily flexural. In the latter case it seems more appropriate 
to write with Ribner (ref. 5) 

-Eil- ♦ ‘Jtt |S a It »{(**. t) 

•fj'v 1 > „ 

where is the flexural damping constant due to the plate and/Sfe*, the 
damping due to the energy radiated to the air. As is shown in Appendix IB 
the noise intensity within the fuselage may or may not be related to the 
acoustic energy radiated by the plate and thus may or may not be 

zero . / 



HA.CA OM llt-20 


9 


We still have to specify the space-wise "boundary conditions on the 
plate and we are led, for the sake of simplicity, to either one of two 
limiting cases. In the first case, the forcing function (the random 
pressure field in the "boundary layer) is characterized "by an integral 
scale so large that at a given time, a skin area "between two stiffeners 
(assumed rigid) is very likely to be subjected to a pressure load of the 
same sign, (see Fig. Ia). This allows us to express y and f as functions 
of t only. 



Fie. lb 


&OR X 


In the second case, the integral scale of the forcing function Is very 
small In comparison with the distance between two stiffeners and the 
behavior of the skin is in the average very much as though the supports 
were removed to infinity (see Fig. Ib) . The real case will in general 
be intermediate between these two limiting examples. However, the first 
case seems to apply to boundary layers of excessive thickness: A 

reasonable guess for the average correlation length might be one dis- 
placement thickness £ *; f or to be larger than the spacing between 

stiffeners (of the order of a foot) the boundary layer thickness £ would 
have to be of the order of five feet or more. This unlikely case is treated 
in reference (7). 

On the other hand, the second limiting case (Fig. Ib) would seem 
to provide a reasonable model for boundary layer thickness not exceeding 
one foot. This is the model discussed now. 

a) The Mean Acceleration 


According to our assumption, the average motion 
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is not sensibly affected by the presence of stiffeners. A large number 
of pulses act on the skin at a given time between two consecutive 
stiffeners . The random pulses may be positive or negative an d thus 
there will be a large number of load reversals between supports . Then 
the effect of the boundary conditions can be expected to become small, 
in the average. Consequently one can define a generalized admittance and 
use it in much the same way as is often done in one -dimensional problems.* 
For instance, the mean square plate displacement is given by 


• % (4t, , hipCo) 

Here the mean square of the forcing function f is related to the 
spectrum (P by: 

CO GO OO 

* ( J[ *F( -fc. > o>) M*.* dk^cLoj 

o •m 


l/jgls the generalized admittance, and kg and CO are respectively 
the wave number in the x direction, the wave number in the z direction 
and the frequency. The determination of ifo is easy once it is realized 
that this expression is the square of the Fourier transform of the 
fundamental solution and so can be written by inspection. Thus an 
average solution of ( 7 ) is: 


and 



m " Tt j„Ju, w )Xk, ikUto 



y ( 4k i, to* 


( 8 ) 


(9) 


Equation ( 9 ) gives the mean square response of an unbounded plate to 
a random forcing function. One should notice that the plate will always 
exhibit resonance no matter what value the damping constant fi may have. 
This resonance occurs, not at a given frequency or at a set of discrete 


* See in particular reference (2). 
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frequencies but over the whole frequency spectrum, whenever the follow- 
ing relationship obtains between frequencies and wave numbers:* 




We can visualize the resonance condition as a crest or ridge In the 
wave space (see Fig. 2) which originates at the line = (S/chV"^ and. 
which becomes higher and steeper as the wave number and the frequency 
increase . Thus the effective damping is a function of the exciting 
frequency. 


b) The Length Scale A 


Equation (6) shows that A 2 is needed as well as (K/a t)2. 
Now A is a length scale. It was defined, say in the x direction as 


MSS < 

£sft>r 




and could he termed the equivalent length of perfect correlation. 


There are various ways of evaluating the integral scale. Perhaps 
the most convenient one for our purpose Is that (found for instance in 
Ref. 2 Eq. 115) which Is derived from the relationship between correla- 
tion functions and spectral functions. Thus if a stationary random func- 
tion J(t) possesses a correlation function which Is sufficiently well 
behaved. 


\ 




* Here resonance is defined as the maximum of the response curve 1/XOO 
holding^ constant. The locus of maximae holding k constant is given 
by: 

V > fa 

{ v. V* V 


These maximae correspond only for zero damping 
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one can define a spectral density function 

— £ (to) = d-T 

^ *o 

Now for the particular case £0=0, this gives 

!«(0- j;?«) Ac 


Since 


Tv* 


JT 

$00 


we have the result that 


X.T till 

2 - 9lo) 


(io) 


We have already obtained a formal representation for the spectral density 
function of the plate. It is the integrand of Eq.. ( 9 ) so we can write, 
in view of ( 10 ) 


.nr to* 

'-at ) 54 H i -I [ V« Ly-^j * y "J hi ,' e 11> 


and a similar expression for Ay. 

* 


cDCoV- 
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III. THE FORCING FUNCTION 


We suppose that a turbulent boundary layer develops on the skin of 
the airplane (on one side of our plate). The forces which excite the 
plate are the pressure fluctuations experienced by the plate itself. 

We assume that all characteristics of the boundary layer are fixed once 
we have specified the boundary layer thickness & , the free stream 
velocity Uao and density^ . In terms of pressure fluctuations this 
implies that at a fixed point of the "wetted" surface of the skin, we 
have for the mean square pressure fluctuations 

and the integral scales, i.e. 





are proportional to S . Also the relative contribution to press ur e 
intensity of the various frequency bands must be a function of Va3,&, 
and only so thatif is a random function of three independent 
variables, x,z,t, F is related to a three -dimensional spectrum by: 




such that: 


TT(<o,-&., k.) - f,u£,fe il 
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Loosely speaking, this means that a characteristic frequency for pressure 
fluctuations is proportional to /J and a characteristic wave length is 
proportional tod . Now the forcing function of Eq_. (7) is a force/unit 
mass so that according to our similarity hypothesis 


i 


2m 





One can thus define a spectral function associated with the forcing 
function f(x,z,t): 







dJz, dM-i-JUo 


and thus 

f*. (^» , 4»«. 7 co) * fjf S'" £. T*. (k. 7 k*.,JZ3 d 2 ) 

<y** ^ l Uo 


Here F 2 is a function of Kq, Kg, and JL only and these are non-dimensional 
variables : 


Jls «o£/u*, 


K. * 

K- 2 . S 3 


-L* 


In terms of these non-dimensional units, Eq. (8) becomes: 




dM, d.Kz dJi Tv.(k.,«», JO 
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and Eq. (9) becomes 



gut [T ( "Jt* K (_ «„ k», JV)<u. <uc>dA 

Cr%hM4 J-=> J-» 


which can he written 


/ dv* 

ITF" 



f/iu. 

cr* h 1 * 



[%■»]* 


(13) 


Again, F 2 is a function of the integration variable only, so that one 
can write 



Equation (l4) yields the two -non -dimensional parameters upon which the 
plate dynamics depend. The first one, C h /SUco is the product of a 


Mach number, (- 


speed of free stream 


-) and a 


^ speed of propagation of waves in the plate 

thickness ratio ( ^m d ar^J^e x thlelfpes^ T he second one 

N plate thickness 

is a non-dimensional damping parameter which is, alike, a function of 
plate and free stream properties. 


5j JtS/lJoo 


If we treat the equation for the integral scale (Eq. 11) in the 
same way, we notice that no new non-dimensional parameter occurs, so 
that, at most 


A „ 
7 " 



Ch_ 

S (Xoo 1 



( 15 ) 
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We now wish, to investigate the form of the functions H and L in 
Eqs. (l^-) and (15) respectively. First, we make an assumption which 
is not strictly necessary hut which simplifies the manipulation of 
Eq.. (12) . We take the function F(K ( , Kj, A) to "be symmetric in K ( 
and f which leads us to define a new wave number. 

K « \l K, i + Kx 


and to write 


Tx ( Kt 5 Kx ® "F ( ^ > ^"0 


Thus, Eq. (13) becomes 


SV-*\ 

\ *u4 

°JL 4 <l/X 1 

0 J 

rz- 

[* K K) 

at J 

1 O’ v h a ' 

> 
/I J 

lo 

?h& k - j 
» «« 




Now, the damping parameter / / Uo ° is assumed small and under these 

circumstances it can be shown (see Appendix II) that 


r<X> 



KdUc 


,0D lid* 


c 1 ^*- oil 

ma i 

L ^ 

• 6X00 *> 



(16) 


The small difference between these two integrals can easily be evaluated 
for arbitrarily small values of /m/u» even though both integrals are 
unbounded as (b -i-0. This leads us to believe that for low damping the 
main contribution to the inner integral comes from the resonance condition 

S 

* 

Thus, if the spectral function F(Jl,K) is reasonably wide, i.e. ifSF/^K«l 
over a large range of K, Eq. (l6) suggests that we write 



T(k,ji) .T (n / 


(17) 
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The requirement that F be flat in K when compared to l/X(K) is equivalent 
to the requirement that the average correlation distance or integral scale 
for the boundary layer pressure fluctuations be small compared to integral 
scale of the plate deflection. Translated in physical terms the simplifica- 
tion suggested here is prompted by the following remark: If the plate 

has some stiffness, it makes little difference whether the forcing 
function is assumed to be distributed over small distances or made of con- 
centrated loads (see Fig. 3) 



FIGUHE 3 


Thus a satisfactory model for the problem at hand would be the impact of 
rain drops on a metal roof. Equation (l 6) allows 1 us to integrate over K, 
to get 



(18) 
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The expression 



can be evaluated only when F(K,Jl) is known. It may be an increasing 
or a decreasing function of <fU«o/ch • In the absence of data on the 
spectral function F, we will not; attempt to define it. 

The function H defined by Eg.. (l4) can be written: 


( 



I SI}* \ 

| 5 Uco 

1 — v = 
lUo J 

Cte) * 

1 IcttJ 


(19) 


where f„ is an unspecified function related by ( 18 ) to the boundary 
layer pressure spectrum. 

In order to determine p ^ 2 we need to find out, in addition, what 
quantities the integral scale \ depends on. Here we make use of con- 
siderations which are similar to those yielding Eq. ( 16 ) (see Appendix III). 
The result Is that 

i' 



where f 2 is another function related to the boundary layer pressure 
spectrum by IIl(^) . How we are able to write Eq.. ( 6 ) as 



Sc 9.' 0»£ 




( 21 ) 


Here 

Ha 

Expression (21) gives the functional dependence of pressure intensity "inside” 
on boundary layer parameters for a typical case. The only q uant ity, not 
immediately available is ). It is probable that we shall have to 

await experimental data to define its numerical value reliably. 
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IV, SPECIAL CASES 


1. Convected Turbulence 


Two authors (ref, 5 & 6) have recently suggested that the boundary 
layer pressure fluctuations at any point of the fuselage skin are caused 
essentially by the passage over the point of a fixed (i.e. time inde- 
pendent) pattern of pressure disturbances carried downstream at a fixed 
convective velocity. So far, experimental evidence in proof or disproof 
is lacking. However, it is interesting to incorporate this special case 
in the general formulation which has been presented. Both the response 
of the plate and the coupling of the plate with the air at rest must 
then be reconsidered. 

a) The coupling of the plate with air at rest in the case of 
convected turbulence 


If a fixed spatial pressure distribution is carried downstream 
on the surface of the plate, it is easy to show that the (infinite) plate* 
response will be of the same kind, i.e., that it will consist of ripples 
which are randomly distributed in space but which travel through the plate 
at the same convective velocity as the boundary layer disturbance. The 
determination of the pressure field inside the fuselage is not in principle 
different for this case and has been carried out in Appendix IB**. The 
result is that for both subsonic moving ripples (with convection velocity 
(Ji < CLl ) and moderately supersonic ones: 



' ' Ja n.(a +Mi 


( 1 . 10 ) 


where 



For higher supersonic speeds, the function of geometry and Mach number 
appearing as an integral is more complicated. The equation (1.10) above 
has the same form as equation (6). On the other hand there is a sharp 
difference in terms of energy radiated by the plate between the subsonic 
and the supersonic case, since no energy at all is radiated by subsonic 
ripples while the supersonic ones do generate some. One must, then, make 


*Here the presence of transversal bulkheads will change the picture 
because of multiple reflections of the ripple. 

**This problem can also be viewed as a steady (randomly bumpy) wing 
problem from the standpoint of a stationary observer. 
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a distinction "between the results in terms of pressure intensity (the 
quantity of practical interest) and in terms of energy radiation. This 
distinction stems from the fact that (as is pointed out on page 7 ) the 
acoustical field investigated is truly a near field. 

b) The response of the plate 

According to the convective hypothesis, time is not an independent 
variable once the convective velocity U, is fixed. Translated in terms 
of the spectral density TT (to, k-^kg) of the pressure fluctuations, this 
means that TC (G) ,k^,k£)is zero, except whenOJ =r or in non- ’ 

dimensional form, when Si ts (Uy W We rewrite equation (12) for this 
special case. v 


-k., A,*, o') = K&.l 

y Uco L UoeJ 


Here £ is the Dirac delta function of the variable «fl . 

Then the plate response becomes 

f W/n \* . „ 4 . r°° 

.(a fc ) " ££ [ 

rCO -qq 


1 - K^F, (k„k 01 U v ^ 

7 (23) 

>0 Jo‘ U ^ 


where 

f| | ) VCv) s 



Now if we assume as before that F is symmetric in % and Kp and 
substitute ^ 


K = Vk, x + 

0:W‘ «±- 

K\ 
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we finally get 


/ dVn x 2, 

{ db ) 


<F' L fr 1, e.* {£ 


tv \ r 


cos 4 *© Fz(£> / cps ©/ 


(24) 


Here 


and 


with 


and 


B - 

W»(3,«os F ( K| 

*■% Ki' K l 

SU “ B^e 


k “ cF" 

The length scaled is indicated "by the following dimensional argument. 

The mean correlation length or integral scale is a weighted average 
of al 1 wave lengths, so that dimensionally 


£ 


/ 


Since resonance dominates the plate response, is given from the 

plate response equation (equation 23) "by the resonance condition 


(uy 

- «■ ■ (tf • 

or 



x -t lu 


(25) 


A similar reasoning would have yielded, in the non -convective case, 
instead of eq. (20). 

Combining (2k) and (25) according to (I. 10) we notice that we can 
still write as in equation (2l). 

fccu <XI <T* h*fb (26) 

Here is a weak function of the Mach number as seen from (1.10). 
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2. The case of zero scale 

Under some circumstances it is possible that the space average of 
the plate motion vanishes, i.e.: 



This does not mean that the normal accelerations at two neighboring points 
show no correlation, but that the 
correlation function becomes negative 
as indicated in Fig. (4) and in such 
a way that its space integral vanishes. 

We can then consider the normal 
accelerations as dipoles rather than 
sources and we are led to a slightly 
different radiation problem. Appendix 
IC shows, however, that if one defines 
a length such that 

x-; _ v? 


The results are again identical in form with those of equation ( 6) . Here 
can be viewed as the mean moment arm of deflection moments. Alternatively 
one can redefine the integral scale as 

Nr e L"i Wf 



where c is a constant. Equation (27) can thus be used to define the 
integral scale in any event. 


( 27 ) 



2k 


NA.CA TM 1^20 


V. SUMMARY OF RESULTS AND DISCUSSION 


Appendix I discusses in addition to the cases mentioned in the text 
a few examples which provide different limiting conditions. Thus the 
observer is brought close to the plate (V«X ). A short time scale is 

considered etc The common feature of all these analyses is that the 

resulting mean noise intensity can always he represented } say hy equation 
(26). We shall therefore retain this equation: 


Ik* (V) £ £* l&<5 

J.a: al a-' 



as the most general statement that we can make at the present time. 
Here 


L 4 

h 

ft 

<x\ 

9o 

c 

U* 

6 

ah 

I* 


% 

Mi 

u, 


= mean square noise intensity inside 
= air density inside 
= speed of sound inside 
= air density in the free stream 
= plate density 
= free stream velocity 
» boundary layer thickness 
= plate thickness 

= viscous damping constant (of units l/time) 

= perpendicular distance between observer and fuselage 
= geometry of the plate 
* Mach number 

= convective velocity of turbulence pattern 
= characteristic velocity in the plate = 




C 
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For all but high supersonic velocities, the dependence, of ^ on is 
quite small and can be disregarded. The function £ (Y,g) , a quantity 
which does not depend on the dynamics of the problem but only on its 
geometry should be modified to take into account the fact that the 
fuselage is a cylinder and not a large flat plate. 

The form of the function S cannot be given here both because no 
information is yet available on boundary layer pressure spectra and 
because S de pends too critically on the type of model assumed. How- 
ever, if [oji is measured while any one of the four variables defining 
S (S , Uac? } c or h) is varied, then the functional form of the noise 
intensity inside a fuselage can be determined. Thus the main contribution 
of the analysis is to diminish the extent of the testing required. 

One of the conclusions which can be drawn from the foregoing equation 
is that unless the boundary layer pressure spectrum is a very sharp 
function of frequency (which would make S very sensitive to ailf/ch ) it 
is not practical to decrease cabin noise by boundary layer suction: Since 

the noise intensity is a weak function of boundary layer thickness, 
decreasing appreciably cabin noise would involve the removal of a 
prohibitive amount of air. 

Another conclusion is that increasing the damping is a very effec- 
tive way of limiting the production of noise of all frequencies, since 
the structure transmits sounds essentially by resonance. 

The analysis which has been presented deliberately omitted some of 
the features of the problem which would influence the results and intro- 
duce new parameters. For instance, the fuselage of commercial airplanes 
is usually subjected to an axial tension as well as other loads. In 
addition the skin Is curved. To account for these features of the prob- 
lem one would introduce further terms In the differential equation 
describing the plate and one could treat It in much the same way as 
has been done here. 

The general methods which have been used are adaptable In addition 
to the study of a germane problem, the fatigue of panels which are 
buffeted by a turbulent boundary layer. 

A NOTE ON TESTING 


The discussion of the various limiting solutions makes it clear 
that for the transmission of boundary layer noise through a structure, 
the ratio of outside (boundary layer) noise to inside (cabin) noise is 
in general a function of boundary layer as well as structural character- 
istics. 

This is to say, first, that an attenuation coefficient cannot be 
defined by testing the structure alone with a standard noise source. 

Thus accurate testing requires at the outset that the model be tested 
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for transmission of a noise similar to boundary layer noise. The main 
property of such a noise ; as we have seen is that it must be random in 
space as well as in time, which precludes the use of one or a few 
concentrated sources as noise generators. The only proper substitutes 
for boundary layer pressure fluctuations are forcing functions whose 
effects on a fuselage are local.* The impact of water drops for instance 
might be found adequate simulation. Further, similarity in testing 
requires the matching of parameters which are ratios of plate and forcing 
function properties. For instance if the forcing function used in the 
test is a turbulent boundary layer, similarity parameters are: 

a*. £ . . 

/ Y ' TUoT ' 



*This is not true of jet noise which is generated away from the fuselage. 



•J 


MCA TM 14-20 


27 


VI. REFERENCES 


1* Lord Rayleigh: The Theory of Sound- Dover Publications, 

New York. Volume II (19^5). 


2. Liepmann, H. W. : Aspects of the Turbulence Problem (Part I) 

ZAMP, Volume III (1952) pp 321-3^2* 


3. Liepmann, H.W. : On the Application of Statistical Concepts 

to the Buffeting Problem. Journal of the Aeronautical Sciences 
19 (1952) pp 793-800. 


4. Timoshenko S. and Young, D. A.: Advanced Dynamics. McGraw-Hill, 

New York (1948). 


5. Ribner, H. S. : Boundary -Layer - Induced Noise in the Interior of 

Aircraft. UTIA Report No. 37 (April 1956) . 


6. Kraichnan, R. H. : Noise Transmission from Boundary Layer Pressure 

Fluctuations. To be published. 


7- Corcos, G-. M. and Liepmann, H- W. : On the Transmission Through 

a Fuselage Wall of Boundary Layer Noise. Douglas Report No. 

sM-19570. (1955) 



28 


NA.CA TM 1420 


APPENDIX I 

TEE RANDOM RADIATION OP A PLANE SURFACE : 
A. FOUR LIMITING CASES 


In order to determine the coupling "between fuselage vibrations and 
cabin air one has to choose a model for the correlation 'iff between the 
normal accelerations at two different points of the plate. The model 
which was discussed and for which equation (4) was made plausible is 
predicated upon two conditions: 


A, That the observer is distant enough so that a large number 
of plate elements vibrating independently contribute sound 
in comparable amounts, i.e. 

\C<Y 

Here as before, ^ is the integral (length) scale for the plate 
normal accelerations and Y is the perpendicular distance between 
the observer and the plate. 

B. That the time scale of the phenomenon is large enough so that 
the differences in phase (introduced by the- unequal distance 
from the point at Y to the various points of a plate element 
of length ?Sm) are unimportant, i.e. 


X « Ai ** 4 


aj_ is the speed of sound in the fuselage air, and 0 is the 
integral (time) scale for the phenomenon: 



Then one can 


choose a simple model for the correlation function 


t 


f ( i *(n) 


( 1 . 1 ) 


where C is the delta function. The normal accelerations are assumed 
perfectly correlated within a length and not at all for distances 
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greater than A • Then 


w 


XT'S? 


4-TT 1 * 


/bVotV -5 

\5tj 


h fix K ’ 


and upon integrating 


_ A 1 - f / 


[f 

+TT 1 * * 

, 2>fc ) . 

U *«■ 


'Which is equation ( 4 ). This case,^ij^, «+ Y 

corresponds to the following conditions. The passenger (or the microphone) 
is far from the plate (in terms of 7^ ), the boundary layer is thick and 
the airplane velocities low. One may well wonder about cases for which 
these conditions do not apply. While it appears difficult to answer such 
a query with generality it is possible to consider other limiting cases. 


For instance let us assume that condition 2 still applies but that 
our observer is extremely close to the plate. This would correspond to 
the following physical case: A thin fuselage skin, a thick boundary 

layer, a low airplane velocity and we are measuring noise by placing a 
microphone very close to the skin and insulating it on all sides except 
the side which faces the skin. Then * 7^« • 



FIGURE $ 

Under these conditions the noise at the microphone is contributed primarily 
from a single plate element which in the average vibrates in phase. The 
evaluation of this contribution is particularly simple, tte can write, very 
nearly 
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If, for the sake of definiteness, we assume the element circular, then 

W” t'iW(^-vr 

and since * Y it is permissible to write 

The pressure intensity is therefore given as 



Thus E q. (6) applies for the very close as well as for the very far 
field when phase effects are not important (T^^ifra^). 

Now assume that we carry on the same e:xperiment but that the boundary 
layer is thin and that the velocity of the airplane is high so that the 
exciting frequencies are high. Let us assume in addition that the skin 
is thick, so that << *y5 0.1 . 



FIGURE 6 
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Nov "the time scale of* the plate motion is short and phase effects are 
prevalent. We define a simple time history in analogy to the space 
description of E<£. (i.l) 



The microphone still receives signals effectively only from one plate 
element and all points within that element vibrate in phase but the 
pressure pulses originating from that element do not arrive at the 
microphone in the same time. Then: 


k- $i($r i j 


-CxN ds.As*. 
/ r, c-u 


( 1 . 3 ) 


A* is simply X, . Equation ( 1 . 3 ) is evaluated by noticing that: 


(i - : 

‘ IldOl 

Here fi are the real roots of g(^) = 0 which are included in the interval 
between and and b. We only have one root, namely ri = r 2 . If we choose 
to integrate, say, with respect to s 2 first we get (assuming again that 
the element is circular 



(1-5) 


and according to (1.4) the inner integral yields: 
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so that 




( 1 . 6 ) 


The time scale appears explicitly in the answer. For the unbounded plate 
however it is simply proportional to S AW just as the time scale for the 
boundary layer pressure fluctuations. 

Finally we may consider a physical case for which phase effects are 
important and for which the microphone has been placed a large distance 
away from the plate, i.e.: b&i « X Y 



BOt/A/OARY 4.AYER 
tm/cka/jess 


3/f/A/ 


O MfCROPHOA/F 


FIGURE 7 
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Now the contrihution from each sub -element of vibrating plate is still 
in the average independent from that of the next one. However, there are 
in addition cancellations from within one element just as in the previous 
case. This will happen if the boundary layer is thin, the airplane 
velocity is high small) and the observer is far from the fuselage wall. 

In order to evaluate this limiting case we first specify the time 
behavior of the correlation function: we write 











( 1 . 7 ) 


so that 

A A 

and we integrate first with respect to . Using the same techniques as 
in the previous example, we get: 

VTC 

i 1 ’* J | ^ (o ; 8,-&v)hy(o> g.-eQ^u ^ 

AO 

Nov we assume that 


l(f (o, #.-9v) mj&p \ 8,(6.-6v^ XLo) 

Integrating with respect to Q^,, 0^, O ^ successively: 




£1 

4TT* 




I 
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For a circular plate of radius R , this would give 


!l 

[ itU(o) X(o-)a; 

\*±L <U 

<TT 

l ) 

) o 1 

S' ! 

i&<ai X In [ i 

* tl 

iTf 1 

1 2>b L 

* 1 


In general, and defining 


a function of the plate geometry and of the distance Y only, we have 




( 1 . 8 ) 
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APPENDIX I 

B. THE NOISE GENERATED BY SKIN RIPPLES OF FIXED .VELOCITY: 


If the turbulence pattern is frozen, as discussed in section IV-l 
ripples will travel through the (infinite) skin at a fixed convective 
velocity. Then the correlation function ifT must be written differently: 


where Ut is the speed of propagation of the ripple (turbulence 
convective speed) and therefore 



The inner integral is of the form 


which can be written 

7 Hl2 

as in part A. The expression 


cWV)* 


( 1 - 9 ) 



56 


MCA TM 1420 


has either one or two real roots depending as M,<1 or t\>\ respectively. 
For , the only real root is 

r, *» r-j. 




i+ 


and thus the inner integral yields 

i 


+ tAtx’i 


so that: 


(K) * El 


O' fi L / MnO f 

4-ir a - ( / J 


rv ( ( itc MjX,*3 


(I. 10a) 


notice that equation (I. 10a) above tends to equation (4) for low 
convective speeds. 


For Hi > | > (1*9) has, in addition to the root 

root given hy 

a*. 

MM 


another 


It is easy to show that this root exists for all values of In order to 
simplify the integration let us assume slightly supersonic conditions; i.e. 
let us write 


M.s i +* 


where 

£ «.* 

Then I 

and £ 


*» S' - iii. 
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and the inner integral = 


1 . „ _j 

'* t M '*' >w 


) 


so that for the supersonic case: 


hi 2 - & V'fcViSs \ h,»» f 

4 K+** ' M « J A.,lc+H,ic l , J 




A result which is save for a constant coefficient the same as 1. 10a. 


(i.lOfc) 
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APPENDIX I 

C. THE GENERATION OF NOISE BI PLATE DEFLECTIONS OF ZERO SCALE 


If the space average of the correlation function is zero and if 
the plate vibrations are isotropic in x* and z' one can define a new 
length scale as a moment arm: 


V'-. [ 

J o 

fff 0 ~ ^*0 ( Si is a fixed point) 

alternatively one can define a modified integral scale 

A" = j“ UJ> U°)Uf 

Here 

X"= «. X 1 

where c is a constant. 

Then, one can idealize the correlation function as 

X’ . X 1 &X 


provided 


It follows that 


K-I&l — 11 4 (.^Sst ) 5 ' l*-*0*(«-*$ 


and integrating with respect to & % 




Now, unless a or b = 0 
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and thus 



If there is a (time) microscale 

( o) **o 

and 



Here the plate has been assumed circular and R is its diameter. 
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APPENDIX II 


THE SIMPLIFICATION OF THE PLATE RESPONSE INTEGRAL 
(Equation 3.6) 


We consider the approximation equation 

(*00 /'CO 

KdK _ 


^ o 




K cLk 


rs^v K^ - n ]*(^,)’ ,n, 


|/LUp 


The right-hand side is clearly unbounded as the damping constant /3 — ** O 
since its value is explicitly proportional to 1 /{3 (see for instance 
Eq.(i8)). On the other hand the difference be twee n- the left and the 
right-hand integrals is finite for {3=0. To show this ve write 


% = 


:h K‘ 


Slice 

Then the left-hand side becomes for {3 - 0 

oL H 


S_ loo 

chu 


f 




and the right-hand side becomes for p = o 


SUL, 
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4ch<C2 


noo 

Jo 
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ir 


trr ; [ 


Ct-|/ C5*'; e C? + D (?-') 


} 


(III) 


(112) 


(113) 


so that the difference D between expressions (III) and (112) is 


or 


P - 
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sa 


ao 


2chll 

<S LL ao 
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(F9 1 + (?'+Q OF) 


8chn 


(mo 



NACA TM U4-20 


hi 


This expression is finite and. of course independent of /3 so that we can 
conclude that the left-hand integral of (l6) is unb o unde d for /3 - 0. .* 
Further, it is clear that D is a regular function of /3 so that the ratio 
of the left-hand side to the right-hand side of Eq.. (l6) can he made ar- 
bitrarily close to unity, by choosing arbitrarily small /3 . If a correc- 
tion is desired a numerical check indicates that Eq.. (Hit-) gives a good 
approximation to the error made even with moderately large dangl i n g. 
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APPENDIX III 

THE EVALUATION OF THE INTEGRAL SCALE 


Our starting point Is Eq. (ll). In terms of non-dimensional variables 
it becomes 


/$!'*] 

\ “k ¥• U«» S 

■» « 

* - 

00 ji u F* (Ji, K, u. ) Ik* 4jl 

bfc J 

' ~ x 

u 


t. 

•SM 


(mi) 


We, now simplify the denominator by writing successively 



C J2 (k*+ u 

S Uco 



Siioo 


(III2) 


Then we define 

£t> -j'- 

$ u«T 

V SlyS'«> x 

(JUtf *•*")* '•C 


Equation (III2) can now be written 

and if we replace Kg by its value at resonance namely 

* The justification for that step is identical to that advanced in 
Appendix II. 
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We can -write 



* % <T'-U i 


Jl* F*(K., 
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